Multipole ordering in /-electron systems on the basis of a j-j coupling scheme 
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We investigate microscopic aspects of multipole ordering in /-electron systems with emphasis 
on the effect of lattice structure. For the purpose, first we construct /-electron models on three 
kinds of lattices, simple cubic (sc), bcc, and fee, by including /-electron hopping through (//cr) 
bonding in a tight-binding approximation on the basis of a j-j coupling scheme. Then, an effective 
model is derived in the strong-coupling limit for each lattice structure with the use of second-order 
perturbation theory with respect to (//f). By applying mean- field theory to such effective models, 
we find different types of multipole ordered state depending on the lattice structure. For the sc 
lattice, a Fag antiferro-quadrupole transition occurs at a finite temperature and as further lowering 
temperature, we find another transition to a ferromagnetic state. For the bcc lattice, a F2u antiferro- 
octupole ordering occurs first, and then, a ferromagnetic phase transition follows it. Finally, for the 
fee lattice, we find a single phase transition to the longitudinal triple-q V^u octupole ordering. 

PACS numbers: 75.30.Et, 71.10.Fd, 75.40.Cx 



I. INTRODUCTION 

It is one of currently important issues in the research 
field of condensed matter physics to unveil exotic mag- 
netic properties of strongly correlated electron mate- 
rials with active orbital degree of freedom. Among 
those materials, in d-electron systems such as transi- 
tion metal oxides, origin of complex magnetic structure 
has been vigorously discussed based on the concept of 
orbital ordering liiSiSi^ Also in /-electron materials in- 
cluding rare-earth and actinide elements, various kinds 
of magnetic and orbital ordering have been foundi^ It 
is now widely recognized that orbital degree of freedom 
plays a crucial role for the emergence of novel magnetism 
in d- and /-electron systems. 

Here we should note that in /-electron systems, spin 
and orbital are not independent degrees of freedom, since 
they are tightly coupled with each other due to the strong 
spin-orbit interaction. Then, in order to describe such a 
complicated spin-orbital coupled system, we usually rep- 
resent the /-electron state in terms of "multipole" degree 
of freedom, rather than using spin and orbital degrees of 
freedom as in c?-electron systems. Among multipole mo- 
ments, there have been intensive and extensive studies 
on dipole and/or quadrupole ordering in /-electron sys- 
tems. In usual cases, magnetic ordering indicates dipole 
one, which can be detected by neutron diffraction experi- 
ments. Ordinary orbital ordering means quadrupole one, 
which can also be detected experimentally, since it in- 
duces lattice distortions due to the spatial anisotropy in 
charge distribution. 

In addition to dipole and quadrupole ordering, in re- 
cent years, possibility of higher-order multipole ordering, 
i.e., magnetic octupole ordering, has been also discussed 
for Cc,Lai_,Bg2&2iiLiiii2JSand Np02,i^'iSiiS*iLima 
to reconcile experimental observations which seem to 
contradict one another at first glance. Very recently, 
a possibility of octupole ordering has been proposed 
also for SmRu4Pi2i^i*2^ It is noted that in these mate- 



rials, crystalline electric field (CEF) ground states are 
Fg quartets with large degeneracy even under a CEF 
potentiali2i2ii2^ In the Fg ground-state multiplet, oc- 
tupoles exist as independent moments besides dipole and 
quadrupole momentsp2& Then, phenomenological theories 
have been developed under the assumption that octupole 
ordering occurs. Note that direct detection of octupole 
ordering is very difficult, since the octupole moment di- 
rectly couples to neither a magnetic field nor lattice dis- 
tortions. However, those phenomenological theories have 
been successful in explaining several experimental facts 
consistently, e.g., induced quadrupole moments in oc- 
tupole ordered states in Cea^Lai-aB u^i^" and Np02>i^*ii 

As mentioned above, thus far, the study on multipole 
ordering in /-electron systems has been almost limited in 
the phenomenological level, mainly due to the complex- 
ity in the treatment of multipole degree of freedom. It 
might be possible to consider a Heisenberg-like model for 
multipole moments, but the interactions among multi- 
pole moments were determined just phenomenologically. 
It is highly required to proceed to microscopic theory, in 
order to understand the origin of multipole ordering in /- 
electron systems. However, it is very hard and practically 
impossible to study multipole ordering in the model re- 
taining all the /-electron states. Then, it is necessary to 
consider a tractable model which keeps correct /-electron 
symmetry. 

One way for such model construction is to use an LS 
coupling scheme. For instance, the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interactions were estimated in 
DyZn,^^ and in CeBg and CeB2C2^ from microscopic 
models using the LS coupling scheme. However, the 
method based on the LS coupling scheme is complicated 
and seems still hard to be extended. One reason for the 
difficulty is that we cannot apply standard quantum-field 
theoretical techniques in the LS coupling scheme, since 
Wick's theorem does not hold. From this viewpoint, it is 
recommended to use a j-j coupling schemed Since indi- 
vidual /-electron states are first defined, we can include 
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many-body effects in systematic ways using theoretical 
techniques developed for the research of d-electron sys- 
temsMi^L^ 

In this paper, in order to investigate how multipole or- 
dering appears in /-electron systems from a microscopic 
viewpoint, we exploit the j-j coupling scheme. We con- 
struct tight-binding models on three kinds of lattices, 
simple cubic (sc), bcc, and fee, by including Coulomb 
interactions among Fg states. In order to discuss multi- 
pole ordering in these models, we derive an effective mul- 
tipole interaction model in the strong-coupling limit for 
each lattice structure by using the second-order pertur- 
bation theory with respect to /-/ hopping integrals, as 
to estimate the superexchange interaction in d-electron 
systems. Then, within a mean-field approximation, we 
clarify what kind of multipole ordering occurs in the effec- 
tive model: For the sc lattice, a T^g antiferro-quadrupole 
transition occurs, while for the bcc lattice, antiferro- 
octupole ordering appears. For the fee lattice with ge- 
ometrical frustration, we find longitudinal triple-q Fs^ 
octupole ordering. 

The organization of this paper is as follows. In Sec. II, 
we introduce a tight-binding model based on the j-j cou- 
pling scheme including only the Fg states. In Sec. Ill, we 
describe the general prescription to derive an effective 
Hamiltonian from the Fg model. In Sec. IV, we show the 
mean-field results of the effective models on sc, bcc, and 
fee lattices. Finally, in Sec. V, the paper is summarized. 



II. HAMILTONIAN 

When we study theoretically the /-electron proper- 
ties, the LS coupling scheme has been frequently used 
to include the effect of Coulomb interactions, spin-orbit 
coupling, and CEF potential. However, as mentioned 
above, it is not possible to apply standard quantum- 
field theoretical technique in the LS coupling scheme, 
since Wick's theorem does not hold. In order to over- 
come such a difficulty, it has been proposed to construct 
a microscopic model for /-electron systems by exploit- 
ing the j-j coupling schemef2S where we include first the 
spin-orbit coupling so as to define the state labeled by 
the total angular momentum j. For / orbitals with an- 
gular momentum £=3, we immediately obtain an octet 
with j=7/2(=3-Hl/2) and a sextet with j=5/2(=3-l/2), 
which are well separated by the spin-orbit interaction. 
Since the spin-orbital coupling is, at least, in the order of 
0.1 eV for / electrons, it is enough to take into account 
the j=5/2 sextet, when we investigate low-temperature 
properties of /-electron compounds in the j-j coupling 
scheme. 

In order to construct the many-body state, we accom- 
modate / electrons in the j=5/2 sextet by following the 
Hund's rule interactions and CEF potential, as we have 
done for c?-electron systems. It has been found that the 
many-electron state obtained in the j-j coupling scheme 
is continuously changed to the corresponding state in 



the LS coupling scheme, as long as those states in both 
schemes belong to the same symmetry group.??' Namely, 
if we based on the spirit of adiabatic continuation, there 
is no serious difference between the states of the LS and 
j-j coupling schemes. Depending on the problem, we can 
use one of the schemes for /-electron systems. For in- 
stance, if we attempt to explain phenomenologically the 
experimental results of /-electron insulators, it is highly 
recommended to use the LS coupling scheme. On the 
other hand, the j-j coupling scheme is rather appropri- 
ate to develop a microscopic theory for novel magnetism 
and unconventional superconductivity of /-electron sys- 
tems. In the present paper, our purpose is to construct a 
microscopic theory for multipole ordering from the view- 
point of spin-orbital complex. Thus, we exploit the j-j 
coupling scheme throughout this paper. 

As described above, we consider only the states with 
j=5/2. The j=5/2 states are further split into F7 dou- 
blet and Fg quartet due to a cubic CEF. In order to con- 
sider multipole phenomena such as octupole ordering in 
/-electron systems from a microscopic viewpoint, in this 
paper we consider only Fg states by assuming large CEF 
splitting energy between F7 and Fg levels. This simpli- 
fication is motivated by the fact that the possibility of 
exotic octupole ordering has been actively discussed in 
Ccj^Lai-ajBg and Np02 with Fg ground state. 

Here readers may be doubtful of the reality of our as- 
sumption, since the Coulomb interaction among / elec- 
trons is naively thought to be larger than the CEF level 
splitting in any case. However, it should be noted that we 
are now considering the /-electron state in the j-j cou- 
pling scheme, not in the original /-electron state with 
angular momentum £=3. As pointed out in Ref. l29l the 
Hund's rule interaction in the j-j coupling scheme is ef- 
fectively reduced to be 1/49 of the original Hund's rule 
coupling. Namely, even if the original Hund's rule cou- 
pling among / electrons is 1 eV, it is reduced to 200 K 
in the j-j coupling scheme. We note that the CEF level 
splitting in actinide dioxides is considered to be larger 
than 1000 K^SiS^ We also recall that the CEF level split- 
ting in CeBg is as large as 500 Ki^ Thus, we safely con- 
clude that our present assumption is correctly related 
to the realistic situation. Of course, in order to achieve 
quantitative agreement with experimental results, it is 
necessary to include also F7 level, since the magnitude of 
the CEF splitting is always finite, even if it is large com- 
pared with the effective Hund's rule interaction. How- 
ever, we strongly believe that it is possible to grasp mi- 
croscopic origin of multipole ordering in /-electron sys- 
tems on the basis of the Fg model, since this model is 
considered to be connected adiabatically from the realis- 
tic situation. We postpone further effort to develop more 
general theory to include all the j=5/2 sextet in future. 

Concerning the /-electron number, in this paper we 
treat only the case with one / electron in the Fg mul- 
tiplet per site. However, this restriction does not sim- 
ply indicate that we consider only the Ce-based com- 
pound. In the j-j coupling scheme, in order to consider 
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FIG. 1: Electron configurations in the j-j coupling scheme 
for CEF ground states, (a) One electron in the Fs for n=l. 
(b) One electron in the Vg for n=3. (c) One hole in the Fg 
for n=3. (d) One hole in the Fg for n=5. 



/"-electron systems, where n indicates local / electron 
number per site, we accommodate / electrons in the one- 
electron CEF levels due to the balance between Coulomb 
interactions and CEF level splitting energy, just as in the 
case of d-electron systems. Thus, the situation with one / 
electron in the Fg multiplet per site expresses both cases 
with n=l in the Fg-Fr [Fig. Q^a)] and n=3 in the Fy- 
Fg [Fig. nib)] systems, where Fa;-Fj, symbolically denotes 
the situation with F^^ ground and Fj, excited states. Fur- 
thermore, we should note that due to the electron-hole 
symmetry in the Fg subspace, the effective model with 
one / electron in the Fg state is the same for that in the 
case with three electrons in the Fg multiplet. Namely, 
the present model also indicates both cases with n—3 in 
the Fg-F7 [Fig. UIc)] and n=5 in the Fy-Fg [Fig. ^d)] 
systems. 

Before proceeding to the exhibition of the Hamiltonian, 
it is necessary to define /-electron operators in Fg states. 
Since the Fg quartet consists of two Kramers doublets, 
we introduce orbital index r (=a and /3) to distinguish 
the two Kramers doublets, while spin index a (=t a-nd J,) 
is defined to distinguish the two states in each Kramers 
doublet. In the second-quantized form, annihilation op- 
erators for Fg electrons are defined as 



/rat = V576ar5/2 + V^ar-3/2, 

for a-orbital electrons, and 

/r/3T = arl/2, 
frpl = ar-1/2: 



(la) 
(lb) 



(2a) 
(2b) 



for /3-orbital electrons, where Urj^ is the annihilation op- 
erator for an electron with the z-component jz of the 
total angular momentum at site r. 

Now we show the Hamiltonian of Fg electrons. For 
the purpose to consider the effective model later, it is 
convenient to express the Hamiltonian in the form of 



H — Hkin + Wloc 



(3) 



where Tikin denotes the kinetic term of / electrons and 
Hioc indicates the local interaction part for Fg electrons. 
In this paper, the kinetic term of Fg electrons is given by 
exploiting the tight-binding approximation. Then, Tikin 



is expressed as 



'^kii 



y ] ^rcr;T'a' frrofr+fiT'. 



(4) 



where /.i is a vector connecting nearest-neighbor sites and 
t'^^.^,^, is the hopping integral of an electron with (r', a') 
at site r-|-/x to the (r, a) state at r. We note that the hop- 
ping integral t'^^.^,^, depends on orbital, spin, and direc- 
tion fi, due to /-electron symmetry. Then, the form of 
the hopping integral is characteristic of lattice structure. 
The explicit form of the hopping matrix will be shown 
later for each lattice structure. Note also the relation 



--tt 



As for the local /-electron term Tiioc, since we assume 
the large CEF splitting energy between Fy and Fg levels, 
it is enough to consider the Coulomb interaction terms 
among Fg electrons. As easily understood from the intro- 
duction of 'spin' and 'orbital' in the j-j coupling scheme, 
the local /-electron term in the Fg quartet becomes the 
same as that of the two-orbital systems for d electrons. 
In fact, after lengthy algebraic calculations for Racah pa- 
rameters in the j-j coupling schemef22i Hioc is given as 



rr r 



(5) 



where rirra =/rrcr/rrcr and nrr=I](T "ttct- ^hc cou- 
pling constants U, U', J, and J' denote the intra-orbital 
Coulomb, inter-orbital Coulomb, exchange, and pair- 
hopping interactions, respectively. These are expressed 
in terms of Racah parameters, and we obtain the rela- 
tion U=U'+J+J' , which can be understood from the 
rotational invariance in orbital spaces^ Note that for d- 
electron systems, one also has the relation J=J'. When 
the electronic wave-function is real, this relation is eas- 
ily demonstrated from the definition of the Coulomb in- 
tegral. However, in the j-j coupling scheme the wave- 
function is complex, and J is not equal to J' in general. 



III. EFFECTIVE MODEL 

In this section, we describe a method to derive an effec- 
tive Hamiltonian by using the second-order perturbation 
theory with respect to hopping integrals. Here we em- 
phasize that the procedure is essentially the same as to 
estimate superexchange interactions for d-electron sys- 
tems, although the calculations are tedious due to the 
existence of orbital degree of freedom. After that, we 
will apply the standard mean-field theory to the effective 
model to depict the phase diagram including multipole 
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ordered states. We believe that it is meaningful to un- 
derstand the complicated /-electron multipole problem 
by using a simple d-electron-like procedure and approxi- 
mations, both from conceptual and practical viewpoints. 

When the Hamiltonian is written in the form of Eq. 
first we solve the local problem 



(6) 



where En denotes the n-th eigenenergy and is the 
corresponding eigenstate with a label a to distinguish the 
degenerate states. Since we accommodate one electron 
per site, the ground state |<i>g) is expressed as 



\n) = 



fPa{r,T,cr) 



|0), 



(7) 



where |0) is the vacuum state of / electrons, a denotes 
the electron configuration in the ground state with one 
electron per site, and Pa{r,T,cr) takes or 1 depending 
on the configuration a. 

Here we consider the formal perturbation expansion in 
terms of '^kin in order to construct the effective model. 
Within the second order, Hctf is generally written as 



'Hcs = 



IcS — ^ / , 



Eq — E„ 



(8) 

where a and b are labels to distinguish the ground states, 
while u is the label for the degenerate excited states. 

Since we consider the situation with one electron per 
site, the intermediate state due to one-electron hopping 
has a vacant and a double occupied site. The dou- 
ble occupied site has six possible states, composed of 
triplet with energy U'—J, doublet with energy 
U'+J{^U-J'), and Ti singlet with energy U+J'. For 
the mathematical completion, it is necessary to include 
all possible excited states in the intermediate process, 
but the calculation becomes complicated. Then, in this 
paper, in order to grasp the essential point of the Ts 
model by avoiding tedious calculations, we include only 
the lowest-energy Fs triplet among the intermediate 
states. This restriction to the intermediate states is val- 
idated, when J is much larger than the hopping energy 
of / electron. Since the /-electron hopping amplitude is 
considered to be small compared with J, even if we also 
include the hybridization with conduction electrons, this 
approximation is acceptable in /-electron systems. 

Let us explain the prescription to derive the effective 
model in the present case. It is convenient to consider 
exchange processes of electrons between two sites, r and 
r'. Since we consider the situation with one / electron 
per site, the initial state |rsi,r's2) is written as 



|rsi,r's2) = /t,jT,jO), 



(9) 



where si and S2 symbolically denote spin and orbital 
states for both electrons. Then, we move one electron 
from the site r' to r. As mentioned above, the inter- 
mediate /^ states at the site r is restricted only as the 



lowest-energy F5 triplet states. Namely, the intermedi- 
ated states should be expressed as \u, r) with the label u 
to distinguish the triplet states, given by 

l + l>r)=/lt/J^TlO)' (10a) 
|0,r) - ifl.fl,^ + fUflf,,)\0)/V2, (10b) 

In order to obtain the effective model Eq. ©, it is 
enough to evaluate the inner product 



Pu, 



(r, u\rs, rs 



(11) 



This quantity is explicitly given by P+i;q|^| = 1, 
Po;aT0i = ^/ V^, and the other non-zero elements are given 

by Pu;s' s — Pu;ss' and — Pu-^rar'a'- 

Then, by including the processes in which an electron 
at r moves first, we obtain the effective Hamiltonian as 

"^cff = — -^S3,S4;si,S2/rS3/rsi/r's^/r's2J (12) 

(r,r') S1-S4 

where (r, r') denotes the pair of nearest-neighbor sites 
and the generalized exchange interaction / is given by 



rr — r 

-'s3,S4:Si,S2 



E\(f'-'^YP* ,P f'^" 

+ it7:lrK-,s,.yPu-,s,,st's-:]/{u' - J). 



(13) 



In order to investigate the multipole ordering, it is 
more convenient to express the effective Hamiltonian 
Eq. (|12|l in terms of multipole operators. For the pur- 
pose, we introduce sonic notations to describe multipole 
operators as 



it = (±V3f" - 

e± = -(f^±V3f^)/2, 



(14a) 
(14b) 
(14c) 
(14d) 
(14e) 



where cr are the Pauli matrices. By using these notations, 
we define one-particle operators at site r as 



A, 



E 

Tt' (7(7' 



frra^Ta-.r' (7' frr'a' i 



(15) 



where A is a 4 x 4 matrix. The multipole operators in 
the Fg subspace are listed in Table |I| 

With the use of above multipole operators, the effective 
Hamiltonian is finally arranged in the form of 



(16) 



where q is the wave vector and Tiiq denotes quadrupole 
interactions. Ti^unq {n=l or 2) denotes interactions be- 
tween Tiun moments and ones between F4„„ and other 
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TABLE I: Multipole operators in the Fg subspace^ The site label r is suppressed in this Table for simplicity. 
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3gu 


3gv 
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Auly 


iulz 


iu2x 


iu2y 


iu2z 


5ux 


buy 
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5av 


5gz 


multipole operator X^'^ 






ol 




■■'y 






•^y 






rp^U 

y 










pseudospin representation 
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(-ay 




fy a^ 







TABLE IL Coupling constants in the eflective model for the 
sc lattice. The energy unit is (l/8)t?/(f/' - J). 
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be 






"3 


6f) 


5(2) 
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-V3 


12 


3 


-4 


-4 
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octupole moments with symmetry different from r4„. 
Tl2q denotes other dipole and octupole interactions. In 
general, 7i2q includes interactions between T4U1 and T^^i 
moments, but we find that such interactions are not in- 
cluded in the models with hopping integrals only through 
(//cr) bonding on sc, bcc, and fee lattices. The explicit 
form of each multipole interaction sensitively depends on 
the lattice structure, as shown in the next section. 

IV. RESULTS 

Now we can calculate the effective interaction between 
two electrons located along any direction r' — r by us- 
ing Eq. (|13l) . if the hopping integral along this direction 
is determined. The hopping integrals of / electrons are 
evaluated by using the Slater-Koster tablei^ In this sec- 
tion, we consider the nearest-neighbor hopping integrals 
through iff a) bonding for three lattice structures, sc, 
bcc, and fee. Then, we present the effective Hamiltonian 
and its ordered states for each lattice. The structure of 
our effective model is consistent with the general form 
of nearest-neighbor multipole interactions on each lat- 
tice derived by Sakai et al^ We follow the notation in 
Ref. |3^ for convenience. 

A. sc lattice 

The nearest-neighbor hopping integrals through (ffcr) 
bonding for the sc lattice are given by 

= [i-77+]ii, (17a) 

^(0,a,0) ^ [J_ (17b) 

= (17c) 

where a is the lattice constant and ti=3(//cr)/14. 

For the sc lattice, the quadrupole interaction term in 
Eq. Ijlfcif) is given by 

Wlq-ai(02VqO^,q^^.+C.p.), (18) 

where c.p. denotes cyclic permutations and C,y—cos{q^a) 
{i'=x, y, or z). The value of the coupling constant ai is 
given in Table UTI 



Note that O^q transforms to {V^Oj^ - 0^^)/2 and 

{—VSOlq — 02q)/2 under c.p. {x,y,z) {y,z,x) and 
{x,y,z) {z,x,y), respectively. The dipole and oc- 
tupole interactions are given by 

H2q = 66[Ti^,"-qTi^:;(C. + Cy) + c.p.], (19) 

and 

H4u„q = 6i"^[Ji™j'r^-+c-p-] 

+ &^"^[^'-q./ir(C. + Cy) + C.p.] (20) 

+bi''\T!^cJtT(c.^Cy) + c.p.)i 

(n) 

where values of the coupling constants bi and b)^ are 
shown in Table ITU 

Note that the form of the hopping integrals Eqs. H17|) 
are exactly the same as those for the Cg orbitals of d elec- 
trons via (dda-) bondingiS2i2& Thus, the effective Hamil- 
tonian has the same form as in the Cg model considering 
only the lowest-energy intermediate states, '^^ when we in- 
terpret that T and a denote Cg orbital and real spin, re- 
spectively. However, the physical meaning of the present 
model is different from that of the Cg model. In particu- 
lar, the effect of a magnetic field is essentially different. 
The dipole moment which couples to a magnetic field H 
is given by J=(7/6)[J*"i-|-(4/7)J'^"2] ^j. ^Yie Tg model, 
while for the Cg model, real spin cr of c? electrons is sim- 
ply coupled to a magnetic field. In contrast to the Cg 
model, a magnetic field resolves the degeneracy in the 
T space even within a mean-field theory for the present 
model, as we will see later. 

By applying mean-field theory to the effective model, 
we find a T^g antiferro-quadrupole transition at a tem- 
perature T=T3g=3ai/kB- As lowering temperature fur- 
ther, we find a r4„i ferromagnetic transition. This ferro- 
magnetic transition can be regarded as a T^u antiferro- 
octupole transition, since the r4„i ferromagnetic state 
with the T^g antiferro-quadrupole moment is equivalent 
to the T^u antiferro-octupole ordered state with the T^g 
antiferro-quadrupole moment. The ground state energy 

is (-3/2)ai - 2be + b[^^ + 2b^^^ per site. 

In Fig. Ofa), we depict an H-T phase diagram. We 
note that the ferromagnetic transition at zero magnetic 
field turns to be a crossover under the finite magnetic 
field. The crossover is drawn by dashed curve, deter- 
mined by the peak position in the magnetic susceptibil- 
ity. Since it is found that the crossover curve is almost 
isotropic in the region shown here, we depict only the 
curve for H || [001]. Note also that under a magnetic 
field, r4„i moments become finite, and then, the T^u 
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2 4 6 8 10 

FIG. 2: Phase diagram and magnetization for the sc lat- 
tice. The Lande g-factor is g,j=Q)/7. (a) H-T phase diagram 
for three magnetic field directions. Solid symbols denote the 
Fsg quadrupole transition. The diamond represents the fer- 
romagnetic transition point. The dashed curve represents the 
crossover to the ferromagnetic state. The definition of the 
crossover is given in the main text, (b) Magnetization as a 
function of magnetic field. 



antiferro-octupole interaction (&6>0) effectively becomes 
a T^g antiferro-quadrupole interaction. Thus, the T^g 
antiferro-quadrupole transition temperature increases as 
H is increased at a low magnetic field region. This be- 
havior reminds us of the experimental results for CeBg, 
although the order parameter in the quadrupole ordered 
phase of CeBg is the T^g quadrupole moment. Magne- 
tization as a function of H is shown in Fig. |2Ib). The 
magnetization is isotropic as H—^O since the r4tji mo- 
ment is isotropic, while anisotropy develops under a high 
magnetic field. 

In Figs.|2ta)-(c), we show specific heat, magnetization, 
and magnetic susceptibility, respectively, as functions of 
temperature. We observe two-step jump of specific heat 
at the quadrupole and ferromagnetic transition temper- 
atures, since we have applied the mean-field theory to 
these second-order transitions. Note that the magneti- 
zation starts to develop below the ferromagnetic transi- 
tion temperature. The magnetic susceptibility exhibits a 



T 1 r 




0.2 0.4 0.6 0.8 1 1.2 



FIG. 3: Temperature dependence of physical quantities in 
the absence of magnetic field for the sc lattice, (a) Specific 
heat, (b) Magnetization, (c) Magnetic susceptibility. 



bend at T^g, while it diverges at the ferromagnetic tran- 
sition temperature. Under the magnetic field, this diver- 
gence turns to be a peak, which defines the crossover to 
the ferromagnetic state in the H-T phase diagram. 



Without magnetic field, the orbital (r) state is con- 
tinuously degenerate in the mean-field theory, although 
such continuous symmetry is absent in this model. As 
has been discussed for an eg electron model such as 
perovskite manganites,'^* quantum fluctuations can re- 
solve this continuous degeneracy, but in the present 
model with the strong spin-orbit interaction, magnetic 
field can resolve this degeneracy. The ground states 
are ferromagnetic with (J^"^) || H, where (• • • ) denotes 
the expectation value. Accompanied O2 ordering is G- 
type [q = (1/2,1/2,1/2) in units of 27r/a] or C-type 
[q = (1/2, 1/2,0)] for H |j [001], while for H || [110], it is 
C-type. For H || [111], there appear C-type O2 ordering 
or equivalent ones in the cubic symmetry. 
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TABLE III: Coupling constants in the effective model for the 
bcc lattice. The energy unit is (882/1327)ti/(f/' - J). 

as a, h h br b['^ 6« 6^ bi'^ bf^ fof bf> 

1 -2 9 2 2 -1 2 -2^3 2 2 -2a/3 



B. bcc lattice 

The hopping integrals for the bcc lattice are given by 

^{a/2,a/2,a/2) ^ _^ + + CT^)/\/3]i2, (21a) 

^(-a/2,a/2,a/2) ^ _^ - - CT")/\/3]i2, (21b) 

^(a/2,-a/2,a/2) ^ _^ + - a^)/\/3]i2, (21c) 

^(a/2,a/2.-a/2) ^ [J _^ - + a")/\/3]i2, (21d) 

where a is the lattice constant and t2—2{f f(T)/21. 

After some algebraic calculations, we obtain the 
quadrupole interaction term for the bcc lattice as 

Wlq = asiOxy-qOxy.q + C.p.)CxCyCz ^^^^ 
+a4[Oyz-qOzx,qSxSyCz + C.p.]. 

The dipole and octupole interactions are given by 

'^2q ^b^xyz^—q^xyZj^^x^y^z 

+h (T5"-qTfq + C.p. )c,c,c, (23) 
+&7[7^l,-q7^,'>s.c,+c.p.], 

and 

T-/ — 74«n jiun i „ ^„ „ „ 

'l-4unq — f^l V^2 — q^2:q ' ^-P- J'^x^y^z 

~^^2 \'^x-qJyq SxSyCz + C.p.] 

Txyz — qi^Jzq ^x^y^z ~l~ C.p.) 

+&i"^[ri^-q(-4r^-«-c, + j,t"s,s.c.) + c.p.], 

(24) 

where Ci/=cos(gya/2) and Sy—sa\(cii,al1\ The values of 

fn) 

the coupling constants a^, 6^, and 6^ are shown in Ta- 
ble EH 

In the mean-field approximation, wc find a r2ti 
antiferro-octupole transition at T2u—'2bz/k^ with 
q=(l,0, 0), and a T^ui ferromagnetic transition at 
a lower temperature. The ground state has the Tc^g 
antiferro-quadrupole moment with the same ordering 
wave- vector as the T2u moment. The ground state 
energy is —as — 65 -f 6^^^ per site. 

In Fig. EJa), we show an H-T phase diagram. Again 
the ferromagnetic transition becomes a crossover under 
the finite magnetic field. The crossover curve determined 
by the peak in the magnetic susceptibility is found to be 
almost isotropic in the region shown here. Then, we show 
only the curve for H || [001]. In the region for high H and 
low T , we find two uniform phases. One is a phase with 



uniform {Txyzr) depending on temperature and another 
is a phase with uniform {Txyzr) which does not depend 
on temperature, as shown in Figs.0Ja) and (c). 

In Fig.QJb), we show magnetization as a function of H . 
We note that the magnetization is isotropic as 7?— s-O as in 
the sc lattice, since the order parameter of the ferromag- 
netic transition is the VauI moment. Note also that the 
jump in the magnetization at [g j ^i^H) / {k^T2u)=^ for 
H II [111] indicates the transition to the uniform state. 

Figures|SJa)~(c) show specific heat, magnetization and 
magnetic susceptibility as functions of temperature, re- 
spectively. We observe two jumps in the specific heat at 
the octupole and ferromagnetic transition temperatures. 
The magnetization begins to develop below the ferromag- 
netic transition temperature. The magnetic susceptibil- 
ity has a bend at T2„ and diverges at the ferromagnetic 
transition temperature. Note that the anomaly in the 
magnetic susceptibility at T2U is very weak. In the pure 
magnetic r2M octupole ordered state, there remains de- 
generacy, while in ordinary magnetic states, degeneracy 
is fully resolved. Thus, the nature of the r2u octupole 
phase is similar to that of the quadrupole ordered phases. 
For instance, the anomaly in the magnetic susceptibility 
is weak at the transition temperature, there is no ordered 
magnetic dipole moment, and another phase transition 
occurs at a lower temperature. 

The ground state is continuously degenerate, since the 
Tiui and Fsg moments arc isotropic in this model within 
the q = (1,0,0) structure. We note that this degener- 
acy is due to the symmetry of the model in contrast to 
the sc lattice. By applying a magnetic field, the ground 
states are uniquely determined. The ground states are 
ferromagnetic phases (J^"^) || H with antiferro Oxy or- 
dering for H II [001] , with antiferro Oyz + Ozx ordering for 
H II [110], and with antiferro Oyz + Ozx + Oxy ordering 
for H II [111]. 

As mentioned in Sec. I, quite recently, a possibil- 
ity of octupole ordering in filled skutterudite compound 
SmRu4Pi2 has been suggested experimentallyi2ii2S In the 
filled skutterudite structure, rare-earth ion surrounded 
by pnictogens form the bcc lattice. Moreover, the 
Fg CEF ground state has been reported in the Sm- 
based filled skutterudite Thus, we expect to apply the 
present model to Sm-based filled skutterudites. When 
we compare our result on the bcc lattice with the ex- 
perimental suggestion, octupole ordering actually occurs 
in our model for the bcc lattice, but r2ii octupole or- 
dered state does not seem to explain the experimental 
results. This discrepancy is due to the suppression of Fy 
orbital, since in the filled skutterudites, conduction elec- 
tron has a„ symmetry, which hybridizes with F7 electron. 
In addition, the level splitting between F7 and Fg is con- 
sidered to be rather small in filled skutterudites. Thus, 
for filled-skutterudite materials, we should consider the 
j=5/2 sextet model in the bcc lattice with itinerant F7 
and localized Fg orbitals. We postpone the analysis of 
such a model in future. 



8 



m 

■"J 



8 

6 
4 
2 

■0 





I 


ft 
H 


1 [001] ^ 
1 [110] • 

1 [111] ■□- 








— • 








'h A 






^•jk - 










14— 

ttl 


0.4 

Tl 


0.6 
'^2u 


0.8 1 



1.8 - (b) 



1 1 1 1 1 r 




H I I [001] 



a I 2 3 4 5 6 ? 





^-0.4 

N 

"0.6 
-0.8 

4 



(c) 




8j[iBH/{kB'2,^) = 5.56 — 



Q 0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 4: (a) H-T phase diagram for the bcc lattice for three 
magnetic field directions. Solid symbols denote the F^u oc- 
tupole transition. The diamond represents the ferromagnetic 
transition point. The dashed curve represents the crossover to 
the ferromagnetic state. As for the definition of the crossover, 
SCO the main text. Open rectangles for H || [111] denote tran- 
sitions to uniform phases, as shown in (c). (b) Magnetization 
as a function of magnetic field for the bcc lattice, (c) Ex- 
pectation value of T2u octupole moment (Txyzr) at each of 
sublattice sites r = (0,0,0) and {a/2, a/2, a/2) under high 
magnetic fields along [111] for the bcc lattice. 
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FIG. 5: Temperature dependence of physical quantities in 
the absence of magnetic field for the bcc lattice, (a) Specific 
heat, (b) Magnetization, (c) Magnetic susceptibility. 



fee lattice 



The hopping integrals for the fee lattiee are given by 



^(0,a/2,o/2) 


= [i + (77+ 


- 4^f^a^)/7]t3, 


(25a) 


^(a/2.0,a/2) 


= [i + ir 


- 4v^f«a^)/7]t3, 


(25b) 


^(a/2. a/2,0) 


= [i + r 


- 4V3f^a^)/7]i3, 


(25c) 


^(0, a/2, -a/2) 


= [1 + 


+ 4V3f^a")/7]t3, 


(25d) 


^(-a/2,0,a/2) 


= [i + ir 


+ 4V3f^a^)/7]i3, 


(25e) 


^(a/2, -a/2,0) 


= [i + (f^ 


+ 4V3f^a")/7]t3, 


(25f) 



where a is the lattice constant and t3={ffa)/8. 

Each multipolc interaction term in the effective Hamil- 
tonian for the fee lattice is given by 

Wlq = ai(02 _q02,qCxCj/ + C.p.) 

+03(02 _qOxy,qSxSi/ + c.p.) (26) 

+ c.p.), 



W2q = &8[Ti^,-qTi^,q(c,C, + C,C,) + C.p.; 

-f 69[r|,"_qT5"qS,s, + c.p.] 



(27) 
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TABLE IV: Coupling constants in the effective model for the 
fee lattice. The energy unit is {l/196)tl/{U' - J). 



ai 
12 



0.3 

64^3 



a4 
192 



195 



b[ 



(2) 



&9 

-336 

1,(2) 



bio 
576 

J2) 



(1) 



6} 



(1) 



-196 



h(2) 



!,(2) 







f,(2) 



224^3 4 193 -336 64^3 2^3 112^3 




and 

n 



4itnq 



+4" 



K-VT(cyc.+c.c,)+c.p.: 

['^K-q'A/q •S^Sj^ + C.p.] 
[Txyz-q{Jzq SxSy + C.p.)] 



(28) 



are 



The values of the coupling constants a^, bi and 5^ 
shown in Table Hvl 

As already mentioned in Ref.jl^, it is necessary to ana- 
lyze the effective model carefully for the fee lattice, since 
the model includes geometrical frustration. It is risky to 
apply directly the mean-field approximation to the effec- 
tive model. First we evaluate the correlation function in 
the ground state using an unbiased method such as exact 
diagonalization on the TV-site lattice. Here we set N=8, 
as shown in Fig. EJa). The correlation function of the 
multipole operators is given by 



Xq^-(l/iV)^e''i-(-- 



(29) 



where (• • • ) denotes the expectation value using the 
ground-state wave-function. 

In Fig. 1^1 (b), we show results for the correlation func- 
tions. The interaction between moments (&io) is 
large, but the correlation function of the moment 
is not enhanced, indicating that the frustration effect is 
significant for an Ising-like moment such as r2u- We find 
large values of correlation functions for J^"^; a-^d 
Oxy moments at q=(0,0, 1). However, there is no term 
in the effective model which stabilizes Oxy quadrupole 
order at q=(0, 0, 1). We note that either of T4U2 and T^u 
ordered states can accompany T^g quadrupole moments. 
Thus, the enhancement of Oxy correlation function in- 
dicates an induced quadrupole moment in r4„2 or 
moment ordered states. Namely, the relevant interac- 
tions are &2^'' and bg, which stabilize the J^"^ and T^" 
order, respectively, at q=(0, 0,1). 

Next we study the ordered state by applying mean-field 

(2) 

theory to the simplified model including only 63 and 
bg. Since the coupling constant bg is slightly larger than 
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Ft 

FIG. 6: (a) 8-site cluster (solid spheres) on the fee lattice 
taken in the calculation, (b) Correlation functions for the 
8-site cluster. The unit of the wave vectors is 2% /a. 



(2) 

b^ , Fsu ordered state should has lower energy than T4U2 
ordered state. The interaction bs stabilizes longitudinal 
ordering of the Fs^ moments, i.e., (T^") || q. 

However, we cannot conclude that the ground 
state is the single-q state ((T^"), (T^^"), (T^")) cx 
{0,0,exp[i2TTz/a]), since there is a possibility of multi-q 
structures. For isotropic moments, single-q and multi-q 
structures have the same energy, and thus, anisotropy in 
the moment is important to determine the stable struc- 
ture. Indeed, the F5„ moment has an easy axis along 
[111] in the Fg subspaceii2*ii In this case, we find that 
a triple-q state is most stable among the single-q and 
multi-q states, since it gains interaction energy in all the 
directions. 

In fact, the mean-field ground-state of the simplified 
model is the longitudinal triple-q Fs^ octupole state with 
four sublattices, i.e., 



5u\ 



(T°r") oc exp[i27r2:/a], 
c>c exp[i27rj//a], 
oc exp[i27rz/a]. 



(Tyr 



(30a) 
(30b) 
(30c) 



This state accompanies the triple-q quadrupole mo- 



men 



{Ozx.) CX (T^;'), 

{Oxyr) (X (T^;')- 



(31a) 
(31b) 
(31c) 
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FIG. 7: The triple-q F^u octupole state. The sur- 
face is defined by 'i'—['^^\tp{6,<l),a)f]^^^ in the polar co- 
ordinates, when the 5/ wave-function is represented by 
'i/{r,8,(j),a)—R{r)tp{9,(j),a), where a denotes real spin. 
White-shift of the surface indicates the increase of the weight 
of up-spin state 0, T)P/''^. 



In Fig. 13 we show symmetry of the charge distribu- 
tion with spin density in the triple-q T^^ octupole state. 
Note that this triple-q structure does not have frustra- 
tion even in the fee lattice. The ground state energy is 
—468 per site, and the transition temperature is given by 
kBT5u—4:bs- We also note that this triple-q T^u octupole 
state has been proposed for Np02 phenomenologicallyii^ 

Let us now evaluate physical quantities in the mean- 
field theory. Figures EJa) and (b) show an H-T phase 
diagram and the magnetic field dependence of the mag- 
netization at r=0, respectively. Note that the magneti- 
zation is isotropic as iJ— >0 due to the cubic symmetry. 
The bend for H || [001] and the dip for H || [110] in 
magnetization indicate transitions to the two-sublattice 
structures. There is anomaly in magnetization also for 
H ]| [111] at the transition to the different sublattice 
structure, but it is very weak. 

Under a high magnetic field, sublattice structures 
change, as shown in Fig. ^ For H ]| [001], we obtain 
a two-sublattice structure with 



(r|«) = 0, 

(rj^r") « exp[i2TTz/a] 



(32a) 
(32b) 
(32c) 



For H 
with 



[110], there appears a two-sublattice structure 



^0, 
^0, 



(T^) = 0. 



(33a) 
(33b) 
(33c) 
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FIG. 8: Phase diagram and magnetization for the simplified 
model on the fee lattice, (a) H-T phase diagram. Solid sym- 
bols denote the Fsu octupole transition, while open symbols 
denote transitions between Fsu octupole ordered states with 
different sublattice structures (see Fig. 0. (b) Magnetic field 
dependence of the magnetization. 



Finally, for H ]| [111], we observe 

(T|-)(xsin[27r(y-z)/a], 
(Ty^r") «sin[27r(z-a;)/a], 
(rf") cxsin[27r(a;-?/)/a]. 



(34a) 
(34b) 
(34c) 



Note also that the triple-q state is fragile under H | 
[110]: (T^^)=0 with a four-sublattice structure for 
g J ^m-qH / [k^T^u) ^ 0.11 at T—Q [this phase boundary is 
not shown in Fig. |HIa)] . 

Figures [rUT al and llOf b) show the temperature depen- 
dence of the specific heat and magnetic susceptibility, 
respectively. At T=T5„, there appear the specific heat 
jump and a cusp in the magnetic susceptibility. In con- 
trast to the sc and bee lattices, there occurs single phase 
transition at zero magnetic field in the case of the fee lat- 
tice. Note also that the cusp structure in the magnetic 
susceptibility is rather strong compared with experimen- 
tal results. ^'^■■^^ Such a quantitative disagreement with 
experiments is considered to originate from the suppres- 
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FIG. 9: Magnetic field dependence of Vsu octupole moments 
at T—Q at eacli of four-sublattice sites: r=(0, 0, 0) (open rect- 
angles), r=(0, a/2, a/2) (open circles), r— {a/2,0, a/2) (open 
triangles), and r=(a/2, a/2, 0) (open diamonds) for the fee 
lattice. 
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FIG. 10: Temperature dependence of physical quantities in 
the absence of a magnetic field for the fee lattice, (a) Specific 
heat, (b) Magnetic susceptibility. 



sion of r7 orbital in our model. The analysis of the j=5/2 
sextet model on the fee lattice is one of future problems. 



V. DISCUSSION AND SUMMARY 

We have constructed Fs models with hopping integrals 
through (//cr) bonding based on the j-j coupling scheme. 
In order to study multipole ordering, we have derived an 
effective model by using the second-order perturbation 
theory with respect to /-/ hopping. By applying mean- 
field theory, we find different multipole ordered states 
depending on the lattice structure. For the sc lattice, 
a Fsg antiferro-quadrupole transition occurs at a finite 
temperature. As lowering temperature further, we find 
a ferromagnetic transition. For the bcc lattice, a 
antiferro-octupole ordering occurs first, and a ferromag- 
netic transition follows it. Finally, for the fee lattice, with 
careful analysis, we conclude the appearance of the single 
phase transition to the triple-q F5„ octupole ordering. 

In this paper, we have not taken into account the ef- 
fect of conduction electron. One may complain about this 
point, since it is believed that the hybridization of / elec- 
trons with conduction electron band is important to un- 
derstand the magnetism of /-electron systems. In fact, in 
the traditional prescription, first we derive the Coqblin- 
Schrieffer model from the periodic Anderson model by 
evaluating the c-f exchange interaction within the 
second-order perturbation in terms of the hybridization 
between /- and conduction electrons4£ Then, we derive 
the RKKY interactions again using the second-order per- 
turbation theory with respect to Jcf. 

In general, the RKKY interactions are orbital depen- 
dent and interpreted as multipole interactions. Such or- 
bital dependence originates from that of the hybridiza- 
tion. Note that the hybridization should occur only be- 
tween /- and conduction band with the same symmetry. 
Here we emphasize that the symmetry of /-electron state 
is correctly included in our calculations. Thus, the struc- 
ture in the multipole interactions will not be changed so 
much, even if we consider the effect of hybridization with 
conduction band, as long as we consider correctly the 
symmetry of / electron states. 

Let us show an example to support our belief. Con- 
cerning the octupole ordering in Np02, we have ex- 
tended the present theory by further including the effect 
of p electrons of oxygen anions^ Namely, we have con- 
structed the so-called f-p model, given in the form of 

Ti = Tif + Tip + Tihyb, (35) 

where TCf and Tip denote the local /- and p-electron 
terms, respectively, and Tihyb is the hybridization be- 
tween p- and /-electrons through (pfa) and (p/tt) . Then, 
it has been found that the structure in the multipole in- 
teractions of the effective model derived from the f-p 
model is qualitatively the same as those obtained in the 
Fg model on the fee lattice. In fact, we have found a 
finite parameter region of Fs^ antiferro-octupole phase. 
Namely, the f-p model on the fee lattice has a tendency 
toward Fs^ antiferro-octupole ordering, which has been 
already captured in the simple (//cr) model. This result 
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suggests that the structure in muhipole interactions is 
determined mainly by the symmetry of /-electron state. 
Most of the effect of hybridization can be included by 
changing effectively (//c) in the multipole interactions 
shown in the present paper. 

However, if the itinerant nature of / electrons is in- 
creased due to the large hybridization and metallicity of 
the ground state becomes significant, the present approx- 
imation inevitably loses the validity and the effect of the 
conduction band should be important. In such a case, 
it is necessary to develop a theory on the basis of the 
orbital-degenerate periodic Anderson model in order to 
include the multipole fluctuations. It is one of future 
tasks. 
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